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ABSTRACT 

Height profiles of mean wind speed and wind velocity fluctuation 
spectra are derived from observations made aboard the R/V Acania while 
anchored in Monterey Bay. The profiles and spectra in conjunction with 
related temperature fluctuation data obtained in a parallel study are 
used to calculate the atmospheric boundary layer turbulence parameters: 
friction velocity, U*, drag coefficient, C^, vertical heat flux (w'T' ) 
and z/L. A general disagreement of the results of this study with each 
other and with previously published results is inconclusive due to small 
size of the data base. The best agreement with results of other studies 
was with the calculation of the drag coefficient using velocity fluctua- 
tion data. 
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I. INTRODUCTION 



In recent years the turbulent structure of the atmospheric surface 
layer has become critically important to scientists working on light 
propagation studies. Turbulence within the surface layer is a primary 
factor in determining such phenomena as the intensity of transmitted 
light when it is received at a point. 

Initial experimental efforts to verify turbulence theory predictions 
were conducted over land. Measurements of wind speed, temperature and 
humidity in those investigations were taken from a stable platform with 
auxiliary instrumentation v/ell protected from external weather. Under 
these conditions relatively accurate measurement of wind speed, tempera- 
ture and humidity could be performed under a variety of weather conditions. 

The marine environment on the other hand presents a number of prob- 
lems, both in accessibility and platform stability. Protection of 
equipment and sensors from damaging effects of the ocean environment is 
also a significant problem. These problems must be solved efficiently if 
scientists are to achieve reliable analyses and predictions of surface 
layer phenomena such as momentum transfer and sensible heat transfer. 

It is a purpose of this thesis to evaluate the R/V Acania (Figure 1) 
as an ocean going platform for turbulence and associated profile data 
collection. The effects of ship's motion and ship's position relative 
to wind direction (i.e. position of sensors relative to the true wind), 
and sensor positions on the ship will be considered. 

Computation of such turbulent transfer parameters as the drag co- 
efficient, Cjj, friction velocity, U*, and vertical sensible heat flux 
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Figure 1. The R/V Acania 
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(w'T') have been made using data collected aboard the R/V Acania during 
1974 and early 1975. The results of these computations will be compared 
with similar results obtained by Cardone (1969) and Businger et al . (1971). 
These comparisons will be used as a basis for evaluating the R/V Acania 
as a platform from which to study turbulence over the ocean. 
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II. TURBULENCE THEORY 



A. GENERAL 

Present boundary layer turbulence theory has its basis in the work 
of Monin and Obukhov (1954). They used a similarity approach to define 
a representative length scale, L, for the surface layer of the atmosphere, 

l = . (,) 

Kg w'T' 

Here g is the gravitational acceleration, T is the ambient temperature, 
and K is the von Karman constant. The selection of the Monin-Obukhov 
Length, L, as a stability scaling parameter is based on the assumption 
that friction velocity, U*, and vertical heat flux (w'T' ) are constant in 
the surface layer. Horizontal homogenity within the surface layer is also 
assumed, and density fluctuations due to pressure changes are neglected. 

Monin and Obukhov (1954) applied this scaling length, using dimension- 
al analysis, to the development of a dimensionless function, 4>^(z/L), 
which can be used to represent the mean horizontal wind variation with 
height, dF/dz, in the surface layer. The following expression is the em- 
pirical relationship for the wind shear in their development. 



dU U* 



f„(2/L) 



( 2 ) 



It can be shown that as vertical turbulent heat flux (w'T') decreases to 
0, indicating neutral stability, (j)^(z/L) must approach 1 if Equation (2) 
is to take on its expected form under neutral conditions. 
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Assuming that convective mixing is negligible under neutral conditions 
it follows that for values of (})^(z/L) near 1 (z<<L) dynamic turbulence is 
of primary importance. Thus the absolute magnitude of L becomes an in- 
dicator of the vertical extent to which mechanical turbulence controls 
the turbulent regime. 



B. STABILITY CONSIDERATIONS 



Observational experiments by Businger et al . (1971) yielded a def- 
inite relationship between Richardson Number, Ri , 



/ 



Ri = 



g Oey/9z) 

Q (8U/3z)^ 



(3) 



and the Monin-Obukhov Length, L, where 0^ is the virtual potential temp- 
erature. Figure 2 from Businger et al . (1971) illustrates this relation- 
ship quite well. The following expressions are approximations of the 
relationship between z/L and Ri , proposed by Dyer and Dicks (1970) and 
Webb (1970) for unstable and stable conditions, respectively. 



z/L = Ri 



(4) 



z/L 



Ri 

1 - a Ri 



(5) 



Here a is an empirically derived constant equal to 0.5. 

It is important to note that Businger et al . (1971) observed the 
Richardson Number to approach a critical value of 0.21 as z/L approached 
+ 0 °. This suggests that as stability increases, the flow becomes es- 
sentially non-turbulent in the surface layer and the effect of mechanical 
turbulence becomes negligible. 
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Figure 2. The dependence of 
Richardson Number on stability. 



Figure 3. Comparison of 
dimensionless wind shear 
observations with interpola- 
tion formulae. 
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The effect of hydrostatic stability on surface layer mechanical tur- 
bulence is further illustrated in Figure 3 from Businger et al. (1971). 
Recalling Equation (2) in light of Figure 3, it becomes apparent that the 
function 4>^(z/L) varies relatively little with respect to z/L for unstable 
conditions approaching a value of 0.5. However, the rapid increase of 
(J)^(z/L) as the atmosphere becomes hydrostaticly stable indicates that the 
neutral assumption breaks down rapidly under increasingly stable condi- 
tions. Discussion in the next section will show that the breakdown of 
the neutral assumption becomes important in the development and use of 
profiles of wind speed versus Ln(z), for the calculation of the turbu- 
lence parameter, U*, the friction velocity. 

C. SURFACE LAYER WIND PROFILES 

Monin and Obukhov (1954) and Panofsky et al . (1960) confirmed that 
under neutral conditions the vertical profile of horizontal wind is loga- 
rithmic. Vertical integration of Equation (2) with (j>^(z/L) equal to 1 
yields 

U(z) = ^Ln(z/ZQ) , (6) 

assuming Uq equal to 0. One must realize, however, that both unstable 
and stable conditions would alter the log profile suggested by Equation 
(6). Figure 4, from Panofsky et al . (1960) illustrates this effect. 

Equation (6) leads to the concept of a surface roughness parameter 
Zq. This parameter is important in determining the shear stress at the 
lower boundary of the surface layer. It might be roughly equated to the 
height of vegetation, for example grass, but it is very difficult to de- 
termine accurately. However, Zq can be eliminated from Equation (6) 
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Figure 4. Effects of unstable and stable conditions 
on height profiles of U. 
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by selecting mean winds, IJi and U 2 at 2 levels, and such that 
(z-|< Z 2 ) and (z-j»Zq). Then inserting appropriate values for U-j and 
^ 2 * derived from a given wind profile, into Equation (6) and subtracting 
the two equations will give 

<<U2 - U J (7) 

U* = ^ ^ , 

Ln(z2/z^ ) 

where k is the von Karman constant equal to 0.35. It is important to 
remember that Equation (7) is only applicable under neutral conditions. 

D. CALCULATION OF U* AND vTT FROM FLUCTUATION PARAMETERS 

In the last section U* was related to the vertical profiles of mean 
horizontal wind in the surface layer. It is also possible to determine 
U* using velocity variance spectra derived from surface layer horizontal 
wind fluctuation data. 

In order to do this it is necessary to have some knowledge of tur- 
bulent kinetic energy balance in wave number space. Kolmogorov's second 
hypothesis plays the key role in developing this knowledge. His hypothe- 
sis states that there is a range of energy wave numbers within which 
Kinetic energy is transferred from one wave number to the next without 
any dissipation. Within this range of wave numbers, called the inertial 
sub-range, turbulent kinetic energy dissipation does not occur. 

The following relationships were proposed by Kolmogorov to define 
the spectral wave number region over which his second hypothesis holds: 



S(k) = C^ e k ■ . 


(8) 


$.^(k) = C2C^k ■ . 


(9) 
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Here S(k) and Sj(k) represent the one dimensional spectra of velocity and 
temperature, respectively, and k is the streamwise wave number. The con- 
stants C-| and have been empirically determined to be 0.5 and 0.25 
respectively. Measurements examined in this study define temporal spectra 
S(f). The streamwise wave number, k, is related to the temporal fre- 
quency, f, by 



where U" is the mean wind at the measurement level. This is based on 
Taylor's (1938) "frozen turbulence" hypothesis. The term "frozen" im- 
plies that the turbulence pattern remains unchanged as it sweeps past the 
probe. 

Of interest in this study is, e, the viscous molecular turbulent 
kinetic energy dissipation. Solving Equation (8) for e produces 



Under neutral conditions it has been shown that turbulent kinetic 
energy production is equal to the rate of molec.ular dissipation of tur- 
bulent kinetic energy, e, hence 



k = ^ 



( 10 ) 



U 




( 11 ) 




( 12 ) 
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Combining Equation (2) for a neutral condition, (J>^(z/L) = 1, with 
Equation (12) and solving for U* gives 




( 13 ) 



where k is von Karman's constant. Thus, under neutral conditions, U* 
can be estimated from either mean wind profiles using Equation (7) or 
from fluctuation data using Equation (13). 

It is also possible to determine the vertical sensible heat flux 
(w'T' ) from measurements of wind and temperature fluctuations and ver- 
tical temperature gradient. Corrsin (1951) proposed the following 

2 

relationship for the temperature structure parameter C-j-: 



Here e has been previously defined and N is the dissipation rate of tem- 
perature variance given by 



on the basis of the temperature variance balance, where 0 is the poten- 
tial temperature. Hence, combining Equations (14) and (15) and solving 
for w'T' yields 



= 3.2 N 



(14) 



N = wT^ 90/9Z 



(15) 



w'j' = _I 

^ ' 3.2 90/9Z 



(16) 
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E. THE BULK AERODYNAMIC METHOD 



One of the comparisons used in this study to evaluate estimates of 
vertical heat flux (w'T' ) in the surface layer was based on work done by 
Friehe (1976). Friehe's results feature the use of bulk aerodynamic for- 
mulae to estimate surface layer heat and moisture fluxes. He used a 
great deal of data from "platforms of opportunity" for these estimates. 
The bulk formulae are necessary because data from these platforms usually 
consist of mean wind and temperature measured at only one level and a sea 
surface temperature. 

The bulk aerodynamic expression for the vertical sensible heat flux 
is as follows: 



(T 3 - T^) ^ (17) 

Here is a heat transfer coefficient constant for the given height 
z-| while iJ-j (m/sec) and , (°C) are the mean wind and temperature at 
and T 3 is the surface temperature. Friehe also considered the sta- 
bility influence on by defining it as follows (for 10 meters): 

= 1 .46 X 10 ■ ^ (U A T > 25 m/sec) , 
or 

= 0.91 X 10 ■ ^ (U A T < 25 m/sec) . 

In Friehe's investigation was 10 meters. 

The next step in Friehe's study was to derive an expression for z/L 
using the bulk aerodynamic expressions for the flux terms in L. He 
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derived the following expression for 1/L: 



L 



-1 



Kg Cm AT 0.61 X 10 ^ Cr AQ 

- {- + 2 

Cn T TT rr 



( 18 ) 



In Equation (18) is the drag coefficient at 10 meters equal to 

_3 

1.32 X 10 , is the moisture transfer coefficient equal to 

_3 

1.3 X 10 , and is the heat transfer coefficient discussed earlier 

in this section. Also AQ (gm/m ) is the moisture difference and AT (°C) 
is the temperature difference between the surface and 10 meters, and 
Ty is the surface mean ambient virtual temperature. 

F. THE DRAG COEFFICIENT 



Having determined values for U* it is possible to calculate a momen- 
tum drag coefficient, C^, that is constant at any given height in the 
surface layer. The drag coefficient is related to U* and the mean wind 
at a given level z-^ by the expression 



C 



D1 
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Many studies have been conducted to determine a representive value 
of Cp at 10 meters (e.g. Cardone 1969). Cardone's summary of C.jq ver- 
sus wind speed for the results of several investigations is presented in 
Figure 5. Calculation of C.|q employing U* determined from both mean 
wind profiles and dissipation rates, when plotted against mean wind at 
10 meters, can be used to evaluate wind observations made aboard the 
R/V Acania. 



22 




gOl. X 



23 



Figure 5. Plot of ten meter drag coefficient, C,«, against mean wind speed at ten meters 
from Cardone (1969). ^ 



III. DATA COLLECTION 



A. PLATFORM 

All data was collected aboard the R/V Acania while anchored in about 
30 fathoms of water in Monterey Bay. The ship's anchorage position, 
shown in Figure 6, lies about 4000 yards northeast of Point Pinos on a 
line between Point Pinos and the Holiday Inn on Del Monte Beach. The 
R/V Acania is 48.4 meters long with a beam of 6.5 meters. She displaces 
224 metric tons with a draft of 2.7 meters. Figure 7 shows exact posi- 
tion of the towers which support the instrumentation and it shov'/s length 
of cable runs and position of the processing and recording equipment. 

While anchored, the ship's engines were used as necessary to main- 
tain ship's head into the wind. This led, on occasion, to placing the 
ship at an angle to the seas, thus exaggerating roll and pitch of the 
platform. Also, turning the screws during data collection probably af- 
fected the sea surface temperature sensor which was placed next to the 
ship at the after end of the deckhouse on the port side in the vicinity 
of the port screw. 

B. INSTRUMENTATION AND PROCESSING EQUIPMENT 

Table I lists the data periods used in this study and indicates 
which instrument levels were operable during each period. 

Mean wind measurements were obtained with Thornthwaite Associates 
cup anemometer wind profile register systems model 104 as shown in Fig- 
ure 8. When rotating, the slotted shaft of the anemometer serves as a 
shutter between a light source and a photoelectric cell. Three cups are 

mounted on stainless steel tubes attached to the main shaft at 120° inter- 
vals# One of these units was placed at each of 4 levels as indicated in 
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Figure 6. The R/V Acania anchorage position. 
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Figure 7. Sensor location on board R/V Acania. 



TABLE I 



*PERIOD 

NUMBER 


DATE 


TIME 


WIND 

LEVEL +,a 


TEMPERATURE 
LEVEL +,a 


1 


3/27/74 


1841 


1 


4 


2 


3/27/74 


1902 


1 


4 


3 


3/27/74 


1928 


1 


4 


4 


3/27/74 


2117 


1 


4 


5 


3/28/74 


1902 


2 


4 


6 


3/28/74 


2059 


2 


4 


7 


8/13/74 


1700 


2 


4 


8 


8/13/74 


1742 


2 


4 


9 


8/13/74 


1742 


2 


4 


10 


9/18/74 


1643 


2 


4 


11 


9/19/74 


1735 


2 


4 


12 


9/19/74 


1846 


2 


4 


13 


9/19/74 


1900 


2 


4 


14 


9/19/74 


1929 


2 


4 


15 


9/19/74 


2011 


2 


4 


16 


9/19/74 


2115 


2 


4 


17 


9/19/74 


2158 


2 


4 


18 


9/19/74 


2243 


2 


4 


19 


9/20/74 


0014 


2 


4 


20 


9/20/74 


0056 


2 


4 


21 


9/20/74 


0140 


2 


4 


22 


9/20/74 


0530 


2 


4 


23 


9/20/74 


0552 


2 


4 
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TABLE I (continued) 



*PERI0D 

NUMBER 


DATE 


TIME 


WIND 

LEVEL +,a 


TEMPERATURE 
LEVEL +,a 


24 


11/21/74 


1800 


2 


4 


25 


11/21/74 


1821 


2 


4 


26 


11/21/74 


1842 


2 


4 


27 


11/21/74 


1903 


2 


4 


28 


11/21/74 


2027 


2 


4 


29 


11/21/74 


2048 


2 


4 


30 


11/21/74 


2109 


2 


4 


31 


11/21/74 


2130 


2 


4 


32 


3/27/75 


1502 


2 


4 


33 


3/27/75 


1519 


2 


4 


34 


3/27/75 


1539 


2 


4 


35 


3/27/75 


1600 


2 


4 


36 


3/27/75 


1644 


2 


4 


37 


3/27/75 


1656 


2 


4 


38 


3/27/75 


1956 


3 


4 


39 


3/27/75 


2020 


3 


4 


40 


3/27/75 


2045 


3 


4 


41 


3/27/75 


2106 


3 


4 



* Hereafter period numbers will be referred to in place of date and time. 



Refer to Figure 7 for instrument level locations. 

® 1 = Levels 2, 3, 4 

2 = Levels 1 , 2, 3, 4 

3 = Levels 1,2,3 

4 = Levels 0, 1, 2, 3, 4 
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Figure 8. 



C.W. Thornthwaite Anemometer Cups. 
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Figure 7. Each was connected to the after deckhouse laboratory by an 
individual coaxial cable, as were all sensors mounted on the 2 towers. 

Placement of the after tower could be a problem. Noting from Figure 
7 the position of the after tower relative to the deckhouse structure, 
it is reasonable to assume that turbulence created by the deckhouse might 
affect at least the third sensor level. 

Data logging for the mean system was accomplished using an NPS de- 
veloped micro-processor based data acquisition system. This MIDAS 
(Microprogranmable Integrated Data Acquisition System) utilizes an 
Intel 8008 central processor to control the sampling, averaging, and re- 
cording of mean meteorological data. All software programming is written 
in PL/M to facilitate the writing of a self-documenting program. 

The operator is interfaced with the system via teletype for full 
duplex input/output communications and program control over the sample 
start-time and the number of samples to be averaged before outputting. 

The operator may also alter the present preset sample list by adding or 
deleting various sensors as they come on line or become inoperable. 

Once initiated, the system is fully automated to sample the tailored list 
of sensors every 30 seconds and periodically print output values averaged 
over the selected interval of from one minute to one hour. In this study 
a ten minute averaging period was used. 

Output values are printed out on the teletype in columnized format 
with the time of print as a leader. The teletype has a paper tape punch 
incorporated which may be activated by the operator to produce a data 
copy concurrent with the print-^out. A magnetic cassette tape recorder 
has been integrated into the system as a third data output device. It 



30 



is planned to interface this cassette to an HP 9830 portable computer 
so that profile and gradient flux estimates can be performed on board 
automatically using the BASIC programming capability of the HP 9830. 
Presently, data cards are punched from the paper tape output and processed 
on the IBM 360 system at NPS. 

Velocity fluctuation measurements were performed using TSI model 
1210 probes with sliding support shields and tungsten wires. The shields 
permitted isolation of the sensing area for the determination of the 
undisturbed velocity, v^, before and after each experiment. The plat- 
inum coated tungsten wire was small enough to resolve the viscous 
dissipation scale without any need to apply wire length corrections. 

The wire was fitted with plated ends for isolating the sensing area and 
thereby minimizing flow disturbance. 

Electronics associated with each probe were a TSI Model 1054B 
linearized anemometer and a TSI Model 1056 variable decade module. The 
anemometer had a linear frequency response from DC to 10 KHz and the 
variable decade module operated with a 0-60 ohm range. 

Sensor placement required exceptionally long cable runs, but any 
decrease in system frequency response had little effect in the frequen- 
cy band of interest. 

Temperature fluctuations were measured using similar sensors. 

The only difference was that a platinum wire was used in place of tung- 
sten. The temperature fluctuation system was designed for a resolution 
of 0.001°C at frequencies up to IKHz. 

The temperature fluctuations were measured using a bridge developed 
by personnel at GTE Sylvania, the GTE Sylvania Model 140. The system 
was slightly modified for use aboard the R/V Acania. 
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The baseband portion of this system is basically a balanced wheat- 
stone bridge excited by a 3 KHz signal with a synchronous detector on the 
output. Segments of a very small diameter platinum wire serve as temper- 
ature sensors in opposite arms of the bridge. The resistance temperature 
coefficients result in an output from the bridge which is proportional 
to the temperature difference between the two probes. 

-4 

The sensor wire is 0.5 centimeters long and 2.5 x 10 centimeters 
in diameter. This extremely small mass allows a response to temperature 
variations of up to 1 KHz, while electronic amplification allov/s tempera- 
ture differences as small as 0.004°C to be observed. 

Both wind and temperature fluctuation data are recorded on a Sanborn 
Model 3950 fourteen channel tape recorder. Real time readout on an eight 
channel chart recorder. Brush Model 240, was used to check the quality of 
the signals coming from the sensors. The charts were also used to select 
the periods to be analyzed for this study. 

Hewlett Packard model HP-2850 temperature sensitive quartz crystal 
probes (Figure 9) v/ere used to obtain mean temperature at the sea surface 
and 4 tower levels. The RF signal from the probes and a reference oscil- 
lator were mixed in a Hewlett Packard HP-2801A readout unit to produce 
a beat frequency v/hose signature can be analyzed to within 0.001 degrees 
centigrade per hertz. Each sensor received pre-experiment calibration 
against a platinum resistance wire thermometer in a temperature con- 
trolled circulating water bath. It was checked over the range of expect- 
ed temperatures. The accuracy in achieving a .005 degree centigrade cor- 
rection factor was a constant for each probe. The tower mounted sensors 
were housed in an aspirated shelter shown in Figure 9 to eliminate 
radiation effects. 
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Figure 9. 

A) C.C. Breidert Company Air-X-Houster Shelter. 

B) Dunmore-type Lithium Chloride Sensor. 

C) Quartz Thermometer Probe. 
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It is felt that placement of the sea surface temperature sensor 
near the side of the ship could have caused erroneous surface tempera- 
tures due to radiation effects near the hull and mixing effects from the 
screws which were sometimes used to maintain ship's position during the 
experiments. The recent movement of this sensor to a position forward 
of the bow should make further sea surface temperature data more re- 
liable. 

Mean relative humidity information was obtained using Dunmore-type 
lithium chloride sensors (Figure 9). This sensor was also placed in the 
aspirated shelter. The basic principle of operation of this sensor is 
resistance change in an electrolytic solution generating a reference 
voltage variance proportional to the relative humidity change. Automatic 
temperature compensation in the instruments meet the following specifi- 
cations for relative humidity: 

+3% humidity below 90% and 
+4% relative humidity above 90%. 

Sensor calibration v/as accomplished by a comparative method using a sat- 
urated saline solution in a closed container. 

C. SPECTRUM ANALYSIS AND INTERPRETATION 

Spectrum analysis is the process through which wind fluctuation data 
is converted to a turbulent energy density spectrum. Appropriate 21 min- 
ute segments of data recorded on magnetic tape were transmitted to a 
Federal Scientific Model UA-500-1 Spectrum Analyzer-Averager. Figure 
10 pictures the analyzer while Figure 11 is a block diagram of the 
analysis process. 

Procedures for converting velocity spectral values, obtained with 
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Figure 10. Federal Scientific Model UA-500-1 Spectrum Analyzer-Average 
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Figure 11. Block diagram of spectrum analyzer. 



an analog spectrum analyzer, to engineering units and for obtaining the 
turbulence parameter e from the velocity spectra are described in this 
section. A parallel discussion for temperature spectra appears in 
Hughes (1976). 

1 . Spectral Plot Scaling 

A necessary procedure was to scale the spectral plots to relate 
RMS input voltages to power spectral densities; variance per unit fre- 
quency. To obtain power spectral density levels, corresponding to RMS 
voltage inputs, calibrated scale charts had to be constructed. 

The charts were constructed as follows. Amplitude scaling was 
accomplished by using an externally generated signal whose RMS value was 
determined by an RMS meter. Using a 0-db (.Iv) input and a spectral gain 
of O-db(xl), a signal with amplitude equal to 0.1 volts RMS was produced 
on a screen display as a spike near the selected frequency. This spike 
was then plotted with an X-Y plotter. Successively, the input gain was 
stepped down to attenuate the amplitude of the input signal by lOdb in- 
crements, and plots of the height of each resulting spectral spike 
(ampl itude) vjere then added to the X-Y plot, at different frequencies. 
These plots represented a graduation of RMS input from a minimum of 0.001 
volts RMS to a maximum of 0.1 volts RMS. An example of such a calibra- 
tion plot on the chart appears in Figure 12. These procedures were 

performed regularly during analyses to insure continual calibration of 
the spectrum analyzer and X-Y plotter. 

For purposes of the chart format the RMS voltages were convert- 
ed to LOG 10 units and a graduated scale was constructed so that the 

logarithm of volts RMS could be interpolated from spectral plots. The 

amplitude scale was adjusted for each spectrum as a function of both 
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Figure 12. Calibration plot for spectrum analyzer. 
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input gain and spectral gain. These values were then converted to PSD 
levels for use in calculating e values. 

The equation (based on Federal Scientific specifications) to 
convert volts RMS to PSD is 



PSD level 



(V^/Hz) = 



0.79(cal. level 
Filter Bandwidth 



( 20 ) 



where the filter bandwidth, 6=0.4, is a function of the analysis range 
selected. 

A more detailed explanation of the spectrum analysis process can 
be found in Lund (1975). 

2. Obtaining e From Scaled Spectra 

The turbulence parameter, e, was obtained from the scaled 
variance spectra on the basis of the universal formula. Equation (8)^ for 
the inertial sub-range in wave number space S(k). This expression pre- 
dicts a -5/3 slope for the spectra when plotted in log-log formats. 

Figures 13, 14, and 15 are typical spectra considered in the 
analyses. It is noted that velocity variance spectra have consistent 
-5/3 slopes. 

Assuming -5/3 slopes for the variance spectra, the intercept 
of the -5/3 slope with the 1 Hz frequency line was the spectral density, 
denoted PSD, value used in computing e. The measured PSD value obtained 
from velocity spectra was converted to a spectral density in engineering 
units by the equation 

S(f) = X PSD 

= (cm/sec/volt)^ x volt ^/Hz = (cm/sec)^/Hz (2T) 
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Figure 13. Typical velocity spectrum with distinct -5/3 slope. Plotter calibration points are 
shown as circles with dots in the center. 
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Figure 14. Typical velocity spectrum with distinct -5/3 slope. 
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Figure 15. Velocity spectrum lacking distinct -5/3 slope. 



where C is the hot-wire calibration factor, the determination of which 
is discussed in the next section. 

Since velocity fluctuations were measured at a fixed point in 
the flow, the resultant spectra are realized at a "temporal" frequency, 
f, and correspond to S(f). 

To obtain an e, the temporal , f, and space, k, scales must be 
related in order to use Equation (8). This is accomplished by using 
Taylor's (1938) "frozen turbulence" hypothesis presented in Section II.D. 
of this study. Thus Equation (8) can also be written as 

fS(f) = kS(k) = k (22) 

2 

where S(f) is a spectral density value with units of (Cm/sec) /Hz. 

The final form of Equation (22) relating e to S(f), which is 
measured, is 

e = ( k^^^fS(f) ^^/^ (23) 

^1 

v/here, empirically, C-j = 0.25 and k = 2irf/Lr. 

From Equation (23), with measured f, SCf), and LT, e can be de- 
termined at each level of interest. This was the final form used to 
determine the e values. 

D. HOT WIRE CALIBRATION 

In-situ calibrations of the velocity sensors were accomplished sev- 
eral times during each experiment. In this procedure recordings were 
made of both the cup anemometer indicated wind speed and the correspond- 
ing hot wire voltage output. The sensor wind speed is given by. 
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= a U + b. 



(24) 



where v is the voltage and U is the mean wind speed for any given level. 
The constants, a and b, are calibration curve slope and intercept re- 
spectively, obtained from the in-si tu calibration. 

To convert the power spectral density levels to velocity units re- 
quires a calibration factor given by 

u' = Cv' (25) 



where C is the calibration factor in cm/sec/volt, v, is the voltage fluc- 
tuation and u' is the velocity fluctuation. Differentiating Equation 
(24) yields 






(26) 



which when substituted in Equation (25) gives 




(27) 



Figure 16 is a sample in-situ calibration curve. 
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Figure 16. Hot wire in-situ calibration plot. 



IV. RESULTS 



A. FRICTION VELOCITY RESULTS 

Friction velocity, U*, was calculated in two ways. In the first 
method (method one) profiles of Ln (z) versus mean wind, U, were plotted. 
Several examples of profiles are given in Figures 17, 18, and 19. 

The best-fit profile for each time period was determined subjective- 
ly. When the sensor at a given level was known to be faulty or was 
highly suspect it was disregarded. When two very different possible pro- 
files existed at any one time, levels 1 and 4 were used to determine the 
profile. Only 4 out of 41 periods were treated in such a manner. On 
only a few other occasions was a particular level considered erroneous. 

In a similar study by Smedely (1975) profiles were drawn using a least 
squares fit to the data. 

After the profiles were drawn, 5 and 10 meter wind speeds were read 
from the profiles for each period and these values were used to solve 
Equation (7) for U*. The results of these calculations for each period 

are shown in Table II as U* . 

P 

Friction velocity was also calculated from fluctuation data (method 
two). First, Equation (23) was solved to obtain a value for the viscous 
molecular turbulent kinetic energy dissipation rate, e, at each level. 
Then Equation (13) was solved to yield U*. Ideally U* at all levels 
would be equal for any given period since U* was assumed to be constant 
in the surface layer. In fact, equipment failure or insensitivity, per- 
iods of nonneutral atmospheric conditions and the finite mathematics of 
computers combined to create less than ideal results. 
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Figure 17. Profiles of mean wind, U, versus Ln(z). 
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Figure 18. Profiles of mean wind, U, versus Ln(z). 
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■igure 19. Profiles of mean wind, U, versus Ln(z). 
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TABLE II 



U * p ( m / sec ) 


U *^( m / sec ) 




^ lOe ^ 


U,|Q(m/sec) 


.21 


CO 


0.82 


1 .50 


7.32 


.26 


.33 


0.88 


1 .40 


8.75 


.27 


.44 


0.78 


2.10 


9.69 


.43 


.42 


2.00 


2.00 


9.45 


.11 


.20 


0.53 


1 .70 


4.82 


.22 


.33 


1.40 


3.20 


5.81 


.16 


.14 


0.90 


0.67 


5.40 


.15 


.16 


0.81 


0.90 


5.33 


.15 


.15 


0.85 


0,85 


5.20 


.14 


.20 


0.88 


1.80 


4.71 


.17 


.25 


0.80 


1.70 


6.00 


.14 


.20 


0,68 


1.40 


5.43 


.10 


.25 


0.40 


2.50 


5.02 


.13 


.18 


0.54 


1 .00 


5.61 


.07 


.20 


0.15 


1 .20 


5.78 


00 

o 


.17 


0.22 


1 .00 


5.38 


.12 


.24 


0.30 


1 .20 


6.87 


.10 


.26 


0.38 


2.60 


5.13 


00 

o 


.19 


0.31 


1 .80 


4.52 


.10 


.13 


1 .40 


2.60 


2.53 


.08 


.13 


1.10 


2.90 


2.41 


.10 


.17 


0.54 


1 .60 


4.30 
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TABLE II (continued 



* PERI 0 D 

NUMBER 


U * p ( m / sec ) 


U * e ( ni / sec ) 






U^Q(m/sec) 


23 


.03 


.16 


0.07 


1 .80 


3.73 


24 


.31 


no 

o 


2.80 


1 .10 


5.90 


25 


.25 


.19 


2.00 


1 .10 


5.63 


26 


.23 


.17 


1 .60 


0.89 


5.70 


27 


.32 


.20 


2.60 


1 .00 


6.28 


28 


.27 


.15 


2.40 


0.76 


5.44 


29 


.19 


.14 


2.10 


1.10 


4.23 


30 


.15 


.14 


1 .40 


1 .30 


3.89 


31 


.15 


.15 


1 .20 


1 .20 


4.34 


32 


.34 


.28 


1.50 


1 .00 


8.75 


33 


.33 


.25 


1 .40 


0.83 


8.70 


34 


.39 


.29 


1 .70 


0.93 


9.50 


35 


.42 


.38 


1 .80 


1 .50 


9.97 


36 


.64 


.39 


3.20 


1.20 


11.36 


37 


.55 


.48 


2.10 


1 .60 


11.90 


38 


.26 


.20 


2.20 


1.30 


5.58 


39 


.38 


.22 


2.90 


0.99 


7.00 


40 


.31 


.22 


3.10 


1 .60 


5.57 


41 


.27 


.22 


2.00 


1 .30 


6.08 



* Refer to Table I for dates and times. 
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Thus when all values of U*^for any given period were within 20% of 
each other a linear average was taken to obtain a representative value 

\ 

for U*. When values at any of the 4 levels differed by more than 20% 
adjacent periods were examined when possible in an attempt to determine 
a trend or isolate a bad level or levels. If no comparison was possible 
the two or three values that agreed most closely were averaged. The re- 
sults of these calculations are listed in Table II as U* . 

Figure 20 is a plot of U* versus U* . It shows little agreement 

P ^ 

between the mean-parameter approach and the fluctuation-parameter ap- 
proach for calculating U*. This lack of agreement might have arisen in 
part from the fact that mean data in this study was analyzed subjective- 
ly while the fluctuation data was analyzed electronically. 

In connection with this lack of agreement, an analytical study of 
possible sensor motion at level 4 resulting from platform roll and pitch 
was conducted. The study showed that vertical motion of the platform 
probably had no significant effect on sensor position until roll or pitch 
or a combination of both exceeded 15°. However, under light wind condi- 
tions particularly, the lateral motion of the anemometers could result 
in erroneously high readings. 

B. DRAG COEFFICIENT RESULTS 

The friction velocity results were used to compute the momentum 

drag coefficient, C-jq. Equation (19) was solved for C-|q at the 10 meter 

level using first ,then U* ^ The results of these computations are 

■ P • ^ 

contained in Table II. Figures 21 and 22 show plots of IIiq versus C-|Qp 
and C-|Q^ respectively. The curve from Cardone (1 969) is plotted for 
comparison. 
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Figure 20. Comparison of friction velocity from Equation (13), 
with friction velocity from Equation (7), U*p. 
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Figure 21. Plot of drag coefficient, C^^p, against mean wind speed, U^q, at 10 meters. 
Average values of C,-, correspond to 1 m/sec intervals. 
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There appears to be credible agreement between Cardone (1969) and 

the C, n results from this study. This lends credence to the U* results 
1 Oe 

over the U*p results in this study. 

C. HEAT FLUX ANALYSES 

Atmospheric surface layer heat flux (w'T‘ ) was calculated using two 
different methods. Method one involved the use of fluctuation data to 
determine the heat flux. Method two employed the bulk aerodynamic ap- 
proach as described by Friehe (1976). 

The first method (method 1) required the calculation of the viscous 
molecular turbulent kinetc energy dissipation rate, e, using Equation 
(23). Then using values for Cy obtained from Highes (1976), Equation (16) 
was solved for (w'T-j ). Table III lists the results of these calculations 
as (w'Tj ) . 

Periods 30, 31, and 37 had values of 0 for (w'T^ ) because 30/30 was 

found to be essentially equal to 0. Here 30/3z was calculated by Hughes 

(1976) from mean temperature profiles drawn for each period considered. 

Period 38 lacked temperature fluctuation data necessary for the calcula- 
2 

tion of Cy used in the calculation of e. 

2 

For all periods, (w'T‘) was calculated for each level where Cy val- 
ues were available. Next a representative value for (w'Tj ) in the sur- 
face layer was obtained by a method similar to that used to obtain a 
mean value of U* in section IV. A. of this study. In any given period 
values within 20% of each other were averaged. Adjacent periods were 
inspected for possible trends when necessary. 

The heat flux calculated with method one was used to solve Equation 
(1 ) with Tq= T.|q + 273°C. Then z/L was calculated at each level for 
every period where (w'Tj) was available and the results were plotted 
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TABLE III 



^PERIOD 

NUMBER 


w'T' X 10^ 
(°C-m/sec) 


w'T^ X 10^ 

{°C-m/sec) 




1 


-0.76 


-2.17 


-3.25 


2 


-0.48 


-2.68 


-3.36 


3 


-1 .84 


-3.06 


-3.47 


4 


-2.33 


-3.19 


-3.71 


5 


-5.29 


-0.56 


-1 .27 


6 


-8.77 


-1.35 


-2.56 


7 


1 .10 


a 


a 


8 


0.62 


a 


a 


9 


0.45 


a 


a 


10 


-0.18 


-0.18 


-0.41 


11 


2.14 


0.25 


0.45 


12 


3.39 


0.32 


0.65 


13 


1 .76 


0.27 


0.60 


14 


3.28 


0.33 


0.65 


15 


0.85 


0.33 


0.63 


16 


0.18 


0.18 


0.37 


17 


2.11 


0.26 


0.42 


18 


1 .14 


0.36 


0.76 


19 


1 .96 


0.58 


1.41 


20 


1 .41 


0.37 


1.60 


21 


1 .13 


0.35 


1.60 


22 


6.94 


0.61 


1 .55 
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TABLE III (continued) 



*PERI0D 

NUMBER 


w'T' X 10^ 
(°C-m/sec) 


w'TI X 10^ 
(°C-m/sec) 


T,o(“c) 


23 


1.83 


0.57 


1 .67 


24 


-1 .40 


0.36 


-0.68 


25 


-3.22 


-0.32 


-0.62 


26 


-4.51 


-0.29 


-0.56 


27 


-7.83 


-0.27 


-0.48 


28 


-8.46 


-0.06 


-0.12 


29 


-26.60 


0.02 


-0.03 


30 


0 


0.02 


0.07 


31 


0 


0.02 


0.05 


32 


-21 .90 


0.53 


0.67 


33 


-10.05 


0.21 


0.26 


34 


68.20 


0.55 


0.64 


35 


32.00 


0.58 


0.64 


36 


1 .25 


0.63 


0.61 


37 


0 


0.37 


0.34 


38 


b 


-0.46 


-0.91 


39 


-2.21 


-0.71 


-0.11 


40 


-0.15 


-0.34 


-0.66 


41 


-4.55 


-0.30 


-0.55 



* Refer to Table I for dates and times. 
v/'T '2 not calculated due to lack of accurate surface temperature 
w'T'i not calculated due to lack of temperature fluctuation data 
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against Ri values obtained from Hughes (1976). Figure 23 is a summary 
of these results. There seems to be little agreement with the curve from 
Businger et al . (1971 ) . 

Method two employed the bulk aerodynamic approach as described by 
Friehe (1976). Mean temperature profiles were used to determine 
(T^ - T^q). Table III lists values of (T^ - T^q) for all periods except 
6, 7, and 8 when surface temperature was considered erroneous. Then 
Equation (17) was solved for (w'T^) with selection of the appropriate 
constant, C^, based on the value of LTaT. Results of these calculations 
are shown in Table III as (w'Tp . 

Figure 24 shows a plot of (w'T‘^ ) versus (w'T'^). The lack of agree- 
ment between ( w ' T ' -, ) and (w'T' 2 ) reiterates the general lack of agreement 
between the mean and fluctuation methods using the data available for 
this study. 

Finally, the bulk method of calculating L was used to develop a se- 
cond estimate of z/L. There Equation (18) was solved for 1/Land then z/L 
was developed for all 4 levels for all periods except 6, 7, and 8, when 
accurate temperature data was lacking. Figure 25 is a summary plot of 
z/L versus Ri where L v;as calculated with the bulk method. Comparison of 
results from this method with Businger et al . (1971) results yields no 
better agreement than the previous case although there appears to be a 
definite trend for the bulk method results. 

Mean temperature profiles used to determine T^q for use in the bulk 
method were developed in the same way that mean wind profiles were de- 
veloped as discussed in section IV. A. of this study. The mean temperature 
values from levels 1 through 4 for each time period were plotted against 
Ln(z). The best straight line fit to these points was drawn subjectively. 
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Figure 23. Comparison of z/L from Equation (1) with Ri . 
Average values of z/L correspond to 0.25 Ri intervals. 
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Figure 24. Comparison of sensible heat flux (w'Tl) from 
Equation (15) with sensible heat flux (w'T^) from Equation (17). 
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Figure 25. Comparison of z/L from Equation (18) with Ri. 
Average values of z/L correspond to 0.25 Ri intervals. 



62 



V. RECOMMENDATIONS AND CONCLUSIONS 



A. RECOMMENDATIONS 

Information on sea conditions would help in determining an average 
ship's motion. This information would be especially valuable under low 
wind conditions, because under such conditions sensor movement could cause 
significant sensing errors, particularly with the anemometers. 

B. CONCLUSIONS 

In conclusion, there is insufficient data at the present time to 
accurately evaluate the R/V Acania as a platform from which to study sur- 
face layer turbulence over the ocean. Many improvements have been made 
in the instrumentation and associated equipment used aboard R/V Acania 
since the data examined in this study were obtained. Future data should 
be more accurate and thus may yield better agreement in comparisons with 
other results. 
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